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I. INTRODUCTION 

Ten years ago, Kane and Fisher showed in their influ- 
ential work on the transmission through barriers in one- 
dimensional (ID) conductors, that at T = an arbitrar- 
ily weak barrier leads to perfect reflection at the impurity 
in the presence of repulsive two-body interaction whereas 
not even an arbitrarily large barrier can hinder the elec- 
trons from perfect transmission through the impurity if 
there is attractive two-body interaction^ Employing the 
bosonic representation of a ID fermionic system, the re- 
sults were obtained from a perturbative renormalization 
group (RG) analysis in the impurity strength, treating the 
weak coupling as well as the dual strong coupling regime. 
Arguing that the RG flow cannot cross the marginal fixed 
point line of zero interaction, the authors obtained their 
universal results independent of the impurity strength. 
The approach and result is closely connected to the per- 
turbational RG treatment of a Brownian particle in a 
periodic potential 2 which experiences a phase transition 
from diffusive to localized behavior at some critical cou- 
pling strength, first observed by Schmidt 

In the following years many authors confirmed and 
also disproved the above results and the controversy was 
only resolved four years ago in favor of the initial re- 
sults by finite-size refermionization at the exactly solv- 
able Toulouse-point X The reason for the steady interest is 
the experimental relevance of the results since ID trans- 
port properties are mainly determined by the transmis- 
sion through the barrier and thus show algebraic behavior 
as a function of the temperature involving non-universal 
exponents, in accordance with Luttinger liquid theoryi^ 
Furthermore, the spectral density of states at the im- 
purity site shows algebraic suppression at the chemical 
potential, but the power-law is governed by a larger expo- 
nent than the bulk exponent of a translational invariant 
model. Recently, evidence was found that the conduc- 
tance in GaAs/AlGaAs quantum wires behaves according 
to the theoretical predictions based on the bosonization 
technique. 6 But it was also pointed out that finite size 
effects might become important since the bosonic repre- 
sentation of the ID Fermi gas is only applicable within an 



asymptotically small regime around the Fermi energy^ 

In this work, we want to apply a novel renormaliza- 
tion scheme to the problem, the flow equation approach 
for HamiltoniansA^ The presented treatment turns out 
to be in a way complementary to the RG approach of 
Kane and Fisher as it is valid for small two-body inter- 
action but arbitrary impurity strength. Further, it does 
not rely on the bosonization technique' 1 ' and it can easily 
be extended to more general one-body potentials. The 
fact that the flow equation approach is a novel renor- 
malization scheme involving different approximations can 
also be seen from the analysis of the two-dimensional 
Hubbard model. Whereas functional renormalization 
schemesii rely on the static limit u> = 0, this kind of 
approximation is not present in the flow equation ap- 
proach since the Hamiltonian formulation of the problem 
is preserved. Unavoidable in the flow equation treatment 
is the neglect of higher order coupling terms in the on- 
site interaction U . Nevertheless, the anti-ferromagnetic 
instability of the strong coupling regime is predicted^ 
Similarly, we will - in the following - recover the exact 
coupling-dependent anomalous dimension even though 
our truncation scheme is valid only for small electron- 
electron interaction. 

The paper is organized as follows. In Sec. II, we will 
first treat the Tomonaga-Luttinger (TL) model within 
the flow equation approach, deriving the exact results 
for correlation functions in the position-time space by 
making explicit use of the bosonization formula. We will 
then employ the usual representation of the fermionic 
field operator, where approximations become necessary. 
Keeping track of only the first non-trivial term, we will 
obtain the algebraic power-law behavior for small two- 
body interaction and the exact anomalous dimension is 
recovered. In Sec. Ill, we will add an impurity to the 
system and show within the same truncation scheme as 
in the case of the TL model that the backward scatter- 
ing constant yields the dominant contribution and cither 
tends to zero or infinity depending on (weak) attractive 
or repulsive interaction. Since the asymptotic behavior 
of the flow equations turns out to be independent of the 
initial conditions, the result also holds for the regime of 
intermediate impurity strength which was not explicitly 
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treated by the RG scheme of Kane and Fisher. We also 
evaluate the spectral density at the impurity after ne- 
glecting the irrelevant forward scattering term. In Sec. 
IV. we close with conclusions. 



II. TOMONAGA-LUTTINGER MODEL 

In the following section, we briefly derive the TL model 
starting from a continuous ID interacting electron gas, 
mainly in order to introduce the notation. We then dis- 
cuss the model via flow equations with and without the 
use of the bosonization formula. The TL model is closely 
related to the n-orbital model in ID on which Wegner 
introduced the method of flow equations.- We also want 
to mention the recent work by Heidbrink and Uhrig on 
ID spinless fermions on a lattice with nearest neighbor 
interaction using continuous unitary transformations^ 



A. The model 

We start with the description of one-dimensional, non- 
relativistic, interacting spinless electrons on a ring of 
length L. The continuum version of the Hamiltonian 
reads 

L/2 

H= J dx^{x)(-^ + U(x))i}(x) (1) 



-L/2 

L/2 L/2 



dxdx'^{x)^{x')V{x - x')ip(x')ip(x) 



-L/2 -L/2 

where the fermionic field ip(x) obeys periodic bound- 
ary conditions and the canonical anti-commutation rela- 
tion. V(x) resembles the two-body potential between the 
fermions and we also included a possible one-body poten- 
tial U(x), which shall be zero for this section. Perform- 
ing the Fourier transformation tp(x) = L^ 1 / 2 ^2 k e tkx c k , 



V(x) = L^V 



with k,q = 2nn/L, n G Z 



and focusing on the low-energy properties in the high- 
density limit of the system, the energy dispersion can 
be linearized around the Fermi points ±/cf where kp = 
2ttN/L, N being the number of electrons. With the 
Fermi velocity vp = kp/m the momentum representa- 
tion of the Hamiltonian is then given by 



H=v F J2\k\c{c k + 1 ± L J2 

k k,k'q 



where c£ creates (annihilates) the plane wave with wave 
number k. Further, we neglected - and will neglect 
throughout this work - contributions from the number- 
of-particle operator, since they vanish in the thermody- 
namic limit. 



v q cl^ q c k cl, +q c k < (2) 



The system can be formulated by two anti-commuting 
species of fermions, the left- and right-movers, c% , re- 
lated to the two Fermi-points ±kp} 4 A mathematically 
more rigorous treatment is possible if the two branches of 
the linear dispersion are extended, ranging from k = oo 
to k = — oo. 15 To assure a well-defined ground-state, the 
states with negative energy have to be filled up with Dirac 
fermions, d k . After proper normalization, the creation 
and annihilation operators of the left- and right-moving 
density fluctuations then obey exact bosonic commuta- 
tion relations. These operators are defined as 



R - „-l/2T r R 



1/2 £< 



k-q u k 



, by — n q 1 I 2 ^ c k+ 



q C k 



where q > and n q = Lq/(2n) E N. The relation be- 
tween the auxiliary and physical fermions is given by 



„R — 



c k k > 
d k k < 



„l _ 



c k k < 
d k fc > 



(3) 



The commutation relation thus reads [b q , b^, ] = 5i ; i>S qtq >, 
with i. i' = L, R. 

The interaction term of the Hamiltonian given in Eq. 
(J2J can now be decoupled into left- and right-moving 
fluctuations if we assume the interaction to be long- 
ranged in position space, i.e. v q /vF <C 1 for q ^> q Cl 
where q c -C kp denotes the interaction cutoff This 
suppresses scattering processes from one branch to the 
other, i.e. t/i-processes according to the g-ology model, 
due to v±2k F ~ 0. A somehow counter-intuitive approach 
is to assume a delta-interaction in position space. A RG 
analysis then shows that repulsive gx-processes scale to 
zero, i.e. are marginally irrelevant^ The so-called g%- 
processes, associated with the density-density interaction 
between left- and right- movers, are then renormalized to 

* 17 
92 = 92 - 01 — 

To also express the kinetic term by the newly intro- 
duced bosonic operators, the Kronig relation has to be 
employed, 



^ Hcfc* - c L k 'c L k ) = £ q (bfb« + 6f ftf) (4) 
<?>o 



ft 



where we neglected contributions of the number-of- 
particle operator of the left- and right-movers^ The re- 
sulting TL model resembles the fixed point Hamiltonian 
of ID interacting electron systems^ 



H 



2irvp 



b^b L q ) 



q>0 

^ E «»« W + W + b f b f + b fbf) 
q>a 



(5) 



Note that the kinetic term is renormalized by scat- 
tering processes which conserve the energy, denoted as 
54-processes in the g-ology models Also in the flow 
equation treatment of the n-orbital model, one has to 
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include scattering terms which conserve the number of 
electrons above and below the Fermi points to the diag- 
onal HamiltonianiS 



B. Solution via flow equation 

Since the Hamiltonian in Eq. J^J is bilinear in bosonic 
operators and the left- and right-moving fluctuations are 
only coupled by the the same wave number, diagonaliza- 
tion is straightforward via Bogoljubov transformation. 
But we want to diagonalize the Hamiltonian via flow 
equations, i.e. we perform a continuous unitary trans- 
formation parametrized by the flow parameter £. 

In differential form, the infinitesimal transformations 
read d(H(£) — [r)(£), H(£)], characterized by the anti- 
Hermitian generator 7](£) — —ry{£). The task is to prop- 
erly define the generator such that the interaction term 
vanishes in the limit £ — > oo . Wegner proposed to choose 
i] = [Ho , H) , where Hq denotes the diagonal Hamilto- 
nian which is assumed to be solved. This choice guar- 
antees that the trace of the squares of the off-diagonal 
elements is a monotonically decreasing function of £, i.e. 
d £ Tr[(H - H ) 2 } < 0. 

But the above choice of ij also induces energy-scale sep- 
aration, i.e. for a simple model Wegner showed that the 
energy difference between the eigenstates of Ho, which 
are just being decoupled, is given by £~ 1 / 2 . Even though 
all generators which diagonalize the system for I — > oo are 
essentially equivalent, a systematic decoupling scheme 
becomes important when approximations are involved. 
For a recent review on the flow equation method, see 
Ref. |U 

In order to abbreviate the notation we define b q = 



and &_ r 



for q > 0. Neglecting the shift of the 



ground-state energy, the 
takes the following form: 



-dependent Hamiltonian then 



h(£) = 5>,(*)&& + \ E + K h \) 



H a + V 



q^O 



(6) 



The initial conditions are given by u° = co q (£ = 0) = 
v F \q\{l + v q /2irv F ) and u° q = u q {£ = 0) = \q\v q /2-jr. The 
generator of the infinitesimal transformations is canoni- 
cally defined as 



r)=[H ,H] 



y^uj q u q (b q b- 



(7) 



where we already used the fact that the flow equations 
preserve the symmetry between left- and right-movers, 



i.e. 



u>q(£) = w~q(£) and u q (£) 



j{£) for all L The 



above choice is likely to eliminate the interaction term V 
for £ — > oofSi but also different generators are possible, 2 - 
The commutator [r/, H] yields the following contribu- 



tions: 

[77, H ] = -J2 ^Wbqb- q + bl q b\) (8) 

fo,V] =- + W (9) 

The flow equations d(H = [r/, H] then read 

d e u> q = -4uj q u 2 , d t u q = ~Auj 2 u q . (10) 

Obviously u> 2 — u 2 — const and with u q {£ = 00) — we 
have u) q = cu q (£ = 00) = (uj 2 (£) — u 2 ^)) 1 / 2 . Inserting the 
initial conditions for £ — yields the well-known result 
u q = v F \q\ y/1 + v q /irv F . 

In order to investigate the flow of observables we need 
to know the ^-dependence of r\ q = —oj q u q . From the 
above equations we obtain 



sgn{v q )Q. 



~,2 



2sinh(4w 2 £ + C 9 ) 



sinh(C g ) 



The explicit ^-dependence of the parameters ui q and u q 
is given by 

oj q (£) = cosh{AE q (£))uj q + smh(AE q (£))u° q (11) 
Uq {£) = cosh(4 J B 9 (£))u° + smk{AE q (l))u° q (12) 

where we defined E q (£) = J Q d£'r] q (£'). One can now con- 
vince oneself that these solutions indeed yield the correct 
boundary values for £ — > 00. 

Since [ri(£),r](£')] = we can also calculate the 
unitary operator U that diagonalizes the TL Hamil- 
tonian. Generally it is given by U = U(£ = 00) 

with U(£) = £exp(/p d£'r](£')) where the operator 
C denotes the ^-ordering operator, defined in the 
same way as the time-ordering operator. We obtain 
U = exp(^^ E*(b q b- q — b'_ q b q )) with the familiar 



relation tanh(4i?*) = — v q / (2nv F - 
E q {l = 00). 



) , where we defined 



To determine the flow of the observables we will use the 
representation of the ID fermionic field operator of the 
left- and right- movers, ^Pl/r{x), which involves the oper- 
ators of the bosonic density fluctuation^^ It is given by 
4>l/r{ x ) — F L / R e' l ' L / R ^ x \ where the operator F L / R low- 
ers the number of left- and right-movers by one, respec- 
tively and commutes with the bosonic operators &q , i.e. 
[b { J\ F L/R ] = 0. 23 The phase field is defined as 



4>l/r(x) 



Ev 1/2 

q>0 



(b 



T9 C 



tf Tq e ±tqX ) 



(13) 



where, for convenience, we omitted the ultraviolet con- 
vergence factor. Due to the preserved symmetry between 
left- and right-movers, only the transformation of the 
right-moving field needs to be considered. 
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Since U(£) only consists of bosonic operators and thus 
commutes with F R , we have U%l) R (x)U^ = F R e u<t>R ^ ui 
for all £. Therefore it suffices to consider the flow equa- 
tions for the phase field, i.e. di<j) R (x) = [r),(j> R (x)]. Dur- 
ing the flow of <fi R (x,£) = J2 q >o <^y( x J ^) different wave 
numbers do not mix and we are allowed to limit ourself 
to the flow of 4>*(x,e) = <pf(£)(tf*(x) + (x) with 

4%' a {x) = {b Tq e^ x - b^e***) . (14) 
This yields the following flow equations: 

d t <p* = -2n 9 <p% , derf = -2 VqV >f (15) 

With the initial conditions ip R (£ = 0) = n q l ^ 2 and 
(fig (£ = 0) = the solution is given by 

<pf(i)=n; 1 / 2 cosh(2E q (t)) , 
rf(£) = -n- 1 / 2 smh(2Eg(£)) . 

In the limit £ — > oo we recover the well-known result for 
the transformed field operator 

-^(b-ge-^-blgj**))) , (16) 

with c q = cosh(2i?*), s q = sinh(2£?*), and 
s 2 = (ujg/Cj q — l)/2, see e.g. Ref. Correlation 
functions in position-time space are now easily calcu- 
lated. 



The initial conditions are given by g(£ = 0) = 1 and 

L I R 

ip q (£ = 0) = 0. The above ansatz implies that we only 
keep track of the first non-trivial term generated by the 
flow equations. 

The flow equations depend on the decoupling scheme 
of the higher order terms. If one neglects bilinear terms 
in the bosonic operators without normal ordering them, 
the parameter g is left un-renormalized, i.e. g — 1 for all 
£. The flow equations then read 

d t <p* = -2 Vg <p*; , 

d t <f$ = -2rwp*-2n- 1 / 2 ri 9 . (18) 

An analytic solution is now possible. In the limit £ — > oo 
the final result is given by 

<p*{l = oo) = n- 1/2 Cg , <f$(£ = oo) = -n- 1/2 s q , 

where we abbreviated c q = c q — 1. 

Expanding also the fermionic field up to linear bosonic 

operators, i.e. ip R (x) -> F R (1 + J2 q >o n q 1 <Pq( x ))i 
we see that the flow equation approach yields the 
correct result of Eq. (|16fl up to terms which are 
linear in the bosonic operators, i.e. Uif} R (x)U' ! — > 

FrQ- + T, q>0 n>q 1/2 (c q <l>f(x) - s q <t> q (x))- Note that the 
Bogoljubov-coefficients s q and c q are exactly recovered. 

Consider now the transformation of the fermionic lad- 
der operators of the right-movers cj?. The flow equations 
dgc R = [77, c R ] do not close and we will again truncate 
the series after the terms which are linear in the bosonic 
operators, 



C. Approximations 

It is not surprising that the flow equation approach 
could be applied in a rather straightforward manner be- 
cause we took advantage of the bosonization technique. 
Now we want to calculate the commutator [rj, tp R (x)] 
directly, employing the usual representation il) R (x) — 
L~ 1/2 J2k eikXc k- With the definition of Eq. (0} this 
yields 



[r],if) R (x)} = 



q>0 



U q {x)Mx) 



4(£) = g(£)c R 



(19) 



9>0 



L I R 

with the initial conditions g(£ = 0) = 1 and ip q {£ 



0) = and where we defined 4> L q /R = {b Tq c H ± g-b\ q c H Tq ) . 

Neglecting terms which are bilinear in the bosonic op- 
erators in their non-normal ordered form leads to the 
same differential equations as for the flow parameters of 
the truncated fermionic field operator. The solution is 
thus again given by g = 1 for all £ and 



This 



Further we have [ij , cf> q ^ R (x)} = —2ri q 4> R ^ L (x). 
shows that the flow equations d(.ip R (x) = [rj,ip R (x)] do 
not close but generate an infinite series of operators. We 
thus make the following ansatz for the fermionic field: 



^ R {x, £) = 1> a (x) \g(£) + (vfWOc) 



W0 9 L (^))] 



q>0 



(17) 



<p R (£ = oo) 



rt(£ = oo) = -n q ; 1 / 2 Sq 



So far we do not know what effect the above trun- 
cation scheme has got on physical quantities. For this 
reason we will now calculate the occupation function in 
momentum space within this approximation, i.e. n R ' 1 = 

(FDS\cf(£ = oo)c R (£ = oo)\FDS), where \FDS) de- 
notes the ground-state of the Fermi-Dirac sea. We obtain 
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the following result: 



R,l 



q>0 



g>0 q q>0 q 



n * — ' n„ 

q>0 q q'>0 q 



+ V -±Q(k F -q-k) + Y ^B(k F + q-k) 

— ' Tin — ' Tl„ 

q>Q q q>0 q 



(20) 

with the step function 0(fc) which equals one for k > 
and zero otherwise. 

Considering k = k F + q with q > 0, only the last 
term of the above equation contributes to n^ 1 . It is 
now convenient to work with a finite interaction cutoff 
q c = 2irn c /L and choose the interaction potential as a 
step function in momentum space, i.e. s = s q for q < 
q c and zero otherwise. Like this, we avoid ultraviolet 
divergences and are able to obtain analytic results. With 
S™=i 1/n - ► hin c + C for n c — > oo, where C is Euler's 
constant, we find in the thermodynamic limit L — ► oo, 
N/L = const the following expression: 



R,l 
l k F +q 



y2^e(q-q)^-s 2 ln(q/q c 



q>0 

>(l- 



(q/Qcr )/2 



(21) 



The last limit was taken in the case of weak coupling, 
where we assumed that the logarithm stems from a 
power-law expansion. In this perturbative regime we 
thus recover the well-known algebraic behavior around 
the Fermi point k F . We also find the remarkable result 
that the above truncation scheme yields the exact anoma- 
lous dimension a = 2s 2 . 

For wave numbers below the Fermi point, i.e. for k — 
k F — q with q > 0, the calculations are not so plane and 
we will have to make further approximations. In the fol- 
lowing, we have to discard J2 q g '>o CqC q > /(n g n g ')(@(k F + 
q - k)Q(k F +q'-k) + e{k F -q-q'-k)- G{k F - q - 
k)Q(k F —q'—k)) and the constant term J2 q >o(.£q+ s q)/ n q- 
The neglect of the first terms can be justified since they 
resemble higher order contributions to the logarithmic 
expansion and by only considering the first non-trivial 
term in the expansion of i/}r(x), these contributions can- 
not be accounted for. The second, constant term implies 
that the neglected operators should be included in the 
flow more adequately. Normal ordering before neglect- 
ing these operators would result in a renormalization of 
g which would also give rise to a constant contribution, 
partially canceling the first one. Nevertheless, an ana- 
lytic treatment is then not possible anymore and we do 
not want to extend the present approach. For a detailed 
discussion on the observable flow in the case of a numer- 
ically solvable, but non-trivial model, see Ref. Izi . 



With these approximations, one then obtains for the 
whole regime in the limit of small coupling 



nf\. = ll2-s m {q){\WQ c ri2 



l k F +q 



(22) 



The pre- factor 1/2 associated with the power-law 
behavior is also recovered from simple perturbation 
theory 19 This factor could not be recovered unam- 
biguously from the calculations of the approximate 
occupation number of the n-orbital models 

Finally, we want to check if also dynamical quantities 
can be recovered from the above truncation scheme. For 
this we will calculate the approximate Green function 
defined as 

iG<:\x,t)= (FDS\ipn ] {x = 0,J? = oo) 

xe iH * t ip R (x,l = oc)e- iH ' t \FDS) , (23) 

with H* = H(l = oo). 

In order to evaluate the time dependence of the 
fermionic field of the right-movers we have to work with 
a constant potential in momentum space, i.e. v = v q 
for all q. This yields a linear energy dispersion for the 
fixed point Hamiltonian, i.e. H* = ^cX^o M^g^? with 
the renormalizcd Fermi or charge velocity v c = v F {\ + 
v/(7rv F )) 1 ^ 2 . With the help of the Kronig relation of Eq. 
(0J, the time dependence of the fermionic and bosonic 
ladder operators in the Heisenberg picture at £ = oo is 
given by c£(t) = cfe"*^* and b q (t) = ^e" 1 !^*. 

The Green function of Eq. (|23|l can now be expressed 
as a function of the conformal variables £ L / R = x ± v c t 
and yields 



iG^ie, e) = iG<>°(t R ) (i - E - a - ^ R )Y 

q>0 Hq 
D 2 



+ E 



q>0 



~2 

fi. e -i<m r 

n„ 



U nq 



Me 



(24) 



L^Y.k^Qikp-k). 



with iG^°{i R ) 

An important consistency check of the above 
truncation schemes is given by the fact that 
nf 1 = ^ l2 l2 dxe- lkx iG^{x,t = 0). But working 
with a constant potential in momentum space, i.e. a 
delta-potential in position space, yields the well-known 
ultraviolet divergences. In the following we will therefore 
introduce an ultraviolet cutoff and also label it with q c . 

A physical observable is given by the occupied den- 
sity of states p^(w) which is observed in photo-emission 
experiments. In the above approximation, it is given by 



<< 1 r \ 
Pr M 



2tt 



dte iujt iG<'\x 



0,t) 



(25) 



Neglecting the same terms as we did in the case of the 
calculation of the occupation number, we obtain with 
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lj = v c kp — 6j in the thermodynamic limit 

p<'V) = e(w) [1 + 2s 2 In(w/(« c g c )] /(2t™ c ) . (26) 

In the limit of small coupling we thus obtain the well- 
known algebraic suppression at the renormalized Fermi 
energy, again governed by the exact anomalous dimension 
a = 2s 2 Note that for the static correlation function 
n^' 1 the anomalous dimension is being accounted for by 
either left-movers (fc > kp) or right-movers (k < kp) 
whereas both branches contribute in equal parts in case 
of the dynamic correlation function p^ l {ui). Finally, we 
want to mention that Eqs. 121fl and H2bj) are also ob- 
tained within RPA perturbation theory^ 

III. TOMONAGA-LUTTINGER MODEL WITH 
IMPURITY 



B. Solution via flow equations 

We will now set up flow equations which will pre- 
serve the form of the Hamiltonian of Eq. I|27|) . This 
strategy has already proved to be successful in the con- 
text of the spin-boson models The generator of the in- 
finitesimal transformations shall consist of three parts 
r\ = rf e + r] F + rj B according to the canonical generator 
rj c = [Hq, H], but with generalized parameters which will 
be determined later, 

rj™=J2v?(b q b- q -bl q bl) , V F = J2^(b q ~bl) 

q^O q^O 

r i B = Y, < + i>y L b q ) - h.c) . 

q=£0 



The two examples of the previous section showed that 
the flow equation approach can yield the exact anomalous 
dimension even within a rather crude truncation scheme. 
We will now add an impurity to the system and employ 
the same truncation scheme as above, i.e. we use the 
fcrmionic representation of the field operators to repre- 
sent the impurity and only keep track of the first non- 
trivial term. 

A. The model 

Departing from the Hamiltonian of Eq. QJ, the one- 
body potential shall now be given by U(x) = Ao^(a;). 
Performing the same approximations that led to the TL 
model and employing the same notation as in the previ- 
ous section, we arrive to the following Hamiltonian: 

q^O q^O 

= H + H ee +Hi , (27) 

with ip = tp(x — 0). The impurity term Hi is usually 
split up into two contributions, i.e. a forward scattering 
and a backward scattering part. With ip = ip^ + tpR they 
read 

Hi = \ F a>y R + ti>y L ) + x B ^y L + ^y R ) 

= Hp + Hb 

The bosonic operators b q and the fermionic field oper- 
ators are thus correlated via the following commutation 
relation: 

[bq^L/n] = [b\^L/R] = -@(Tq)riq 1/2 il>L/R 

The forward scattering contribution Hp can be ex- 
pressed by the bosonic density fluctuations, i.e. H F = 
L^ 1 J2 g =io n V 2 ^q(b q + b q ), where we neglected the term 
proportional to the number-of-particle operator of the 
left- and right-movers. The initial conditions are given 
by Af {£ = 0) = A and \ B (£ = 0) = A . 



Notice that rj is not normal ordered in the usual sense 
such that the fermionic creation operators are left from 
the fermionic annihilation operators. For the following, 
we thus want to define normal ordering such that the 
creation (annihilation) operators 6j (b q ) are left (right) 

from the 2fcF-operators ip R ipL and ipj^ipR, respectively. 
By this, we take on the standpoint that the bosonic op- 
erators are independent entities with no internal struc- 
ture. 

To clarify the approach we will now present the several 
contributions of the commutator [77, H\. The commutator 
[r/ ee ,H] yields the following contributions: 



lv ee ,H ] 


= 2j2" q r)?(b g b- q + bl q bl) 


(28) 




= 2j2v?u q (b q bl + blb q ) 


(29) 




= 2L- 1 Y,n\' 2 r ]q «\ F _ q {b q + b\) 

q^O 


(30) 


lv ee ,H B ] 


= 2j2Vq e ^ B sgn(q)n- 1/2 (y q + h.c.) 
-2j2 V - e X B n q 1 ^y L + h.c.) 

q^O 


(31) 



In Eq. (|31[1 we defined the operators ^f q = {b q ip R i'L — 
YR^Lbq), whose contributions to the flow will be sup- 
pressed by a proper choice of r\ B . The relevant contribu- 
tions of [77^ , H] are given by 

[ V F ,H ]=J2v^q(bq + bl) , (32) 

q^O 

[ V F ,H ee ]=J2v q u q (b q + bl) . (33) 

q^O 

The commutator [r/ B , H] yields the following contribu- 
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tions {[if, H F ] = 0): 

[r, B ,H }^-^2r,fu, q (^ q + h.c.) (34) 

[ V B , H ee ] « $>*,u, - + h.c.) (35) 

+ 2^r?f u g sgn(g)n- 1 / 2 (^L + 

[r/ B ,ff B ] = 4^f A fl s gn ( 9 )n- 1 /V^K4^L (36) 
+ 2/>X>? A B sgn(g)n7 1 / 2 (^V'Ji + 

q^O 

+ 2pJ2vg^ B {bl^UR ~ V-IV-l) + h.c) 

q^O 

In Eq. I2U we introduced the ^-dependent parameter 



«#0 " 9 



The symbol in Eqs. ijSH) and 
signifies that we neglect normal ordered operators which 
consist of two bosonic operators and one 2fcp-operator 
(ipRtpL or iPl,iPr)- This is our truncation scheme which 
should be valid for small electron-electron interaction v q 
as was pointed out in the previous section. 

We will now look more closely at Eq. (J2EJ). First, we 
note that the factor p = N/L emerges from the regular- 
ization of the anti-commutation relation {i[>l/r> yl/r} ~ 
L~ zl k ~ N/L, having in mind the restriction of the 
sum to the first Brillouin zone. Secondly, we see that new 
interaction terms are generated which can be expressed 
by the bosonic density fluctuations of the left- and right- 
movers, i.e. H R = L~ x E^oK^JV + R lA b « b «' + 
bq'bq)) with general coupling matrices R q ^ q >- This contri- 
bution had to be included into the flow of the Hamilto- 
nian. However, there is evidence based on investigations 
of the TL model with open boundaries that these terms 
represent boundary effects which - in the thermodynamic 
limit - do not alter the universal features of the models 
The evolution of Hr shall thus be neglected. 

We will now specify the flow. The generators rf e 
and r] F are determined according to the canonical gen- 
erator r\ c = [H ,H], i.e. rf q e = -uj q u q and rj F = 



J q"-q 



L X q . The generator rj is chosen in order 
to eliminate the newly generated terms that contain the 
operator (ty q + h.c). This yields 



rif = S ga(q)n-^ 2 X B - 



2^ e 



(37) 



The constant k, which stems from normal ordering, acts 
as a low-momentum cutoff and can be neglected in the 
asymptotic regime since k — ► £~ x / 2 . For £ = 0, we obtain 
k oc vq c where v is a measure of the electron-electron 
interaction and q c denotes the interaction cutoff. For 
attractive interaction, v < 0, and large |u| and q c , this 
might lead to a singular expression of Eq. (|37|l . For 
the following discussion, this parameter regime has to be 
excluded. 



With deH = [rj,H], the flow equations for the bulk 
parameters read 



= -4:U> q u q 



dgu q = 



-4ujgU q 



(38) 



These are the same flow equations as in the case without 
the impurity. The solution is thus given by Eqs. 
and |Q ■ The asymptotic behavior of the coupling u q {£) 
turns out to be crucial for the following analysis. It can 
be deduced from the exact solution and yields u q (£) — > 
u q {u q /Lj q )e- 4 ^ e for I oo. 

The flow equations for the scattering terms read 



de\ B = X B ^ 



4:U)qU q 



8u) q u q 



1 



9>0 

F (, ? 
q - -q K 



dp\ F = -X q (a)* + 3u q u q ) 



+2 P (x B fY: 



q'>0 



UJqi +U q , + K 



(39) 



(40) 



Setting k = 0, we observe after some algebra that the flow 
of the impurity strength delicately depends on the sign 
of the interaction potential v q . For attractive electron- 
electron interaction v q < the absolute value of the 
backward scattering potential X B decreases whereas it 
increases for repulsive interaction. We will now show 
that |A B | tends to either zero or infinity which resembles 
the main result of Kane and Fisher. 

Eq. is easily integrated and yields X B {£) = X B (£ = 

0) exp( ll d£'S{£')), where we defined 



S 



'4uj q U q 



E 

q>0 



Sw q u\ 



-AGiit 



(41) 



and the limit was taken for £ — > oo. If we choose the 
electron-electron interaction potential as a step function 
with momentum cutoff q c , i.e. v q — vQ(q c — \q\)r^ then 
the sum can be performed in the thermodynamic limit, 



4cD^- 

LU 



dqqe 



'), 



where we defined uj = lim 9 ^ u = IhUg— >o u q /q an d 

u) = liniq^o ^q/q- The asymptotic behavior of the back- 
ward scattering impurity potential is thus given by 



X B -> X B £</ 2 



(42) 



where a* B = v/(2itvf + v) and X B denotes the £- 
independent asymptotic part of X B . This leads to the 
scenario described above and we want to stress that the 
asymptotic behavior is independent of the initial con- 
dition, i.e. holds for arbitrary initial impurity strength 
A . In Fig. ^ the qualitative flow of |A B | depending on 
the sign of the two-body interaction v is shown as it fol- 
lows from the one-loop perturbative RG analysis of Kane 
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v<0 v=0 v>0 v<0 v=0 v>0 

FIG. 1: The qualitative flow of |A S | depending on the sign of 
the two-body interaction v as it follows from perturbative RG 
analysis (l.h.s.) and from the flow equation approach (r.h.s.). 
The solid lines resemble analytic results whereas the dashed 
lines are inferred. 



want to mention that a spatially dependent impurity po- 
tential U(x) already couples the forward and backward 
scattering potentials within the above truncation scheme. 

If we set the forward scattering potential X F zero from 
the beginning and simply neglect the normal ordered con- 
tributions of the generated operators (tyg + h.c), the flow 
equations close involving only the generator rj ee . With 
X F = for all I, we obtain the same asymptotic flow 
equations for u> q , u q , and A B and thus the same scenario 
as above. This simplification will be useful for the eval- 
uation of correlation functions, outlined in the next sub- 
section. In this context, we also want to mention Refs. 

where the results of Kane and Fisher are con- 
firmed using Monte Carlo and Bethe ansatz techniques, 
respectively. 



and Fisher (l.h.s.) and from the flow equation approach 
(r.h.s.). The solid lines resemble analytic results whereas 
the dashed lines are inferred. 

Note that the phase separation is due to the asymptotic 
flow in which the system is already almost decoupled, 
i.e. the two-body interaction u q (£) is exponentially small. 
The localization phenomena is thus a genuine low-energy 
phenomena. Under these considerations it appears plau- 
sible that for real (finite) samples this asymptotic regime 
may be hard to reachi 

What is still left to verify is that the forward scattering 
potential will be negligible compared to A B for I — > 
00 . The asymptotic behavior of the inhomogeneous part 
of Eq. is given by 



U q ^ + u ^ 



4ui— 



= ir 1 / 2 -r 1 / 2 etf(2u> qc V£) , (43) 



^ Jo 

, u 
-( 



with the error function erf(x) = 2n 1 I 2 f£ dx'e x — > 

1 — e~ x / (xir) for x — > 00. The differential equation of 
Eq. (|40|l can thus be solved in leading order as £ — > 00 
to yield Af = TT^ 2 p(u/u;)(Xf /Q q ) 2 l a B-V* , 

Since a* B — v/(2ttvf + v), the forward scattering im- 
purity strength tends to zero as £ — ► 00 for v < 2-kvf- 
For v > 2itvf - one must question whether the trunca- 
tion scheme still holds for this coupling regime - A^ goes 
to infinity, but with a weaker power-law behavior than 
the backward scattering potential. The forward scatter- 
ing term will therefore always resemble the less relevant 
contribution to the asymptotic flow. 

Analyzing the above equation for finite system size 
L, we find that the maximum value of |A F | increases 
exponentially as a function of the initial value Ao before 
it eventually tends to zero. Since an extended flow 
equation scheme would couple the forward scattering 
and backward scattering potentials, this behavior could 
result in a significant change of the absolute value of X B 
in finite-size systems which is found in Ref. 0- We also 



C. Spectral density at the impurity 

To discuss correlation functions, the observable has 
to be subjected to the same sequence of infinitesimal 
transformations that led to the diagonalization of the 
Hamiltonian. The flow equations for the observable read 
diipa{x) = [77, i/jr(x)] and differ for x — and x ^ 0. At 
x = 0, the leading terms are given by 



+4'L 



9 + hii>y>L 

q>0 
q>Q 



which is symmetric in the left- and right-movers and only 
the initial condition indicates the evolution of the right- 
moving field. The leading expansion of the field operator 
at x ^ is given by 



1>r(x, I) = ^r(x) \g(x) + J2 

<2>0 



(45) 



The different expressions are crucial in order to distin- 
guish between the bulk and the boundary regime. 

The flow equations for the field operator turn out to 
be rather complicated. This is a general feature of the 
method since flow equations are designed to diagonal- 
ize the Hamiltonian and not to yield simple expressions 
for the observable flow. 24 We will therefore allow for the 
above mentioned simplification of the Hamiltonian flow 
by neglecting the forward scattering contribution to the 
Hamiltonian. Since now only the generator rf e is in- 
volved, we can use the calculations of the previous sec- 
tion to determine correlation functions, but we loose the 
distinction between the bulk and the boundary regime. 
The truncation scheme of the field operator is thus given 
by Eq. I|17fl for au x - Still, we obtain different expres- 
sions for the approximate Green function, as defined in 



9 



Eq. H23|) . f° r attractive and repulsive interaction due to 
the different fixed point Hamiltonians H*. 

For attractive interaction, the fixed point Hamiltonian 
H* is given by the free bosonic bath and the fermionic 
ground-state \FDS) is obtained with the help of the Kro- 
nig relation of Eq. Q). Following the procedure of the 
last section, we thus recover that the spectral function is 
governed by the exact bulk exponent a = 2s 2 throughout 
the system even at x = 0. 

For repulsive interaction, due to the remnant back- 
ward scattering potential, the fixed point Hamiltonian 
has to be transformed further to assure a simple ground- 
state. First, we again map the non-interacting bosonic 
bath onto the non-interacting Fermi gas with linear dis- 
persion employing the Kronig relation. The external po- 
tential with zero width and zero transmission can now 
be absorbed by changing the boundary conditions of the 
physical field operator from periodic to fixed or open 
boundary conditions. To describe a linearized ID sys- 
tem with fixed boundary conditions, only right- movers 
need to be introduced^ We thus transform the original 
ladder-operators of the left- and right-movers according 
to the following substitution scheme: 



-ic k 



IC-k 



(46) 



The sign-change of the momentum in the case of the 
left-movers expresses the fact that there are only right- 
movers; the relative phase factor stems from the condi- 
tion that the physical field operator has to obey fixed 
boundary conditions, i.e. ip(x = 0) = 4>r{x — 0) + 
i> L (x = 0) = 0. 

For finite system size, the transformation is more com- 
plicated since a periodic system of size L can only be 
mapped onto an open system of size 2L in order to ac- 
count for the different wave numbers and to conserve the 
density of the right-movers^ Since we are only inter- 
ested in results valid in the thermodynamic limit, we will 
stick to the simple substitution scheme of Eq. (|46[) and 
divide the transformed fixed point Hamiltonian by two 
in order to preserve the density of the right-movers, i.e. 

Performing the substitution also in Eq. I|17|). one sees 
that only the combination (p q = (ip^ + (fig) enters. With 
the same approximations as in the previous section and 

— 1/2 

with ip = riq (p q for q < q c and zero otherwise, the 
spectral density of states at the impurity is given by 



Q(u>) 



2ttv c 
9(6) 

2nv c \v c q, 



1 + (2^* + <p* 2 ) ]n(Cj/v c q c ) 



with u> = v c kp —uj and v = (<p* + 1) 2 . As usual, for small 
interaction algebraic behavior is inferred and the asterisk 
denotes the limiting value as I — + oo. 

The differential equation for (p follows from Eq. i|18|) 
and yields the final result 



^ = ^ = o)=n- 1 / 2 (e- 2 ^-l) . (47) 

We thus recover the exact boundary exponent v = 
e~ 4E i = l/g, where the last definition corresponds to 
the notation of Kane and Fisher. 



\VrQrJ 



IV. CONCLUSIONS 



In the present work, we applied flow equations for 
Hamiltonians to the TL model with impurity in the 
limit of weak electron-electron interaction but arbitrary 
impurity strength. We formulated the TL model by 
collective density fluctuations which obey Bose statistics, 
but to represent the impurity potential no use of the 
bosonization technique was made. We obtained the 
phase diagram of Kane and Fisher and could provide 
analytic evidence that it also holds in the regime of 
intermediate impurity strength. Simplifying the flow of 
the Hamiltonian by neglecting the forward scattering 
contribution in the initial Hamiltonian, flow equations 
for the field operator could be analyzed to yield the exact 
anomalous dimensions for the different phases as they 
follow from the bosonization approach. Generalizations 
of the external one-body potentials are straightforward. 
The presented approach might also be useful for other 
systems and resembles an alternative and in a way 
complementary method compared to functional renor- 
malization schemes. 
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